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Abstract 

We derive an equation of state (EOS) for strange matter, starting 
from an interquark potential which (i) has asymptotic freedom built 
into it, (ii) shows confinement at zero density (ps = 0) and deconfine- 
ment at high ps, and (in) gives a stable configuration for chargeless, 
/^-stable quark matter. This EOS is then used to calculate the struc- 
ture of Strange Stars, and in particular their mass-radius relation. 
Our present results confirm and reinforce the recent claim ||, ||] that 
the compact objects associated with the x-ray pulsar Her X-l, and 
with the x-ray burster 4U 1820-30 are strange stars. 
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Strange stars (SS) are astrophysical compact objects which are entirely 
made of deconfmed u,d,s quark matter {strange matter). The possible exis- 
tence of SS is a direct consequence of the conjecture p, |J that strange matter 
(SM) may be the absolute ground state of strongly interacting matter rather 
than 56 Fe. Since this hypothesis was formulated, strange stars have been 
studied by many authors ||, ||, 0, f|, but they remained purely theoretical 
entities. This situation changed in the last few years, thanks to the large 
amount of fresh observational data collected by the new generation of x-ray 
and 7-ray satellites. In fact, recent studies |2[] have shown that the com- 
pact objects associated with the x-ray pulsar Her X-l, and with the x-ray 
burster 4U 1820-30, are good strange star candidates. 

Most of the previous calculations of SS properties used an equation of 
state (EOS) for strange matter based on the MIT bag model. In this phe- 
nomenological model the basic features of QCD {i.e. quark confinement and 
asymptotic freedom) are built in. The deconfinement of quarks at high den- 
sity is not obvious in the model. Some people have also expressed more 
general reservations about the bag model and, in particular, about the use 
of a large value for the bag constant ||. In fact, to get SS configurations 
with a mass-radius (MR) relation in agreement with the semiempirical MR 
relations [|T], [T(J for the two SS candidates mentioned above, one must indeed 
use large values of the bag constant B ~ 110 MeV/fm 3 (which in different 
units correspond to B = {hc) 3 B ~ (170.5 MeV) 4 , which is still small com- 
pared with the estimate from QCD sum rules ||, |TT|. But this value is large 
when compared with the "standard value" B = 56 MeV/fm 3 ((144 MeV) 4 ) 
which is able to reproduce the mass spectrum of light hadrons and heavy 
mesons [12]. Large values of B are also not allowed by the requirement SM 



is stable in bulk [|H| [TJ]]: e.g. for m s = 150MeV the permitted values of the 
bag constant are in the range (56 - 78) MeV/fm 3 (see Fig.l of ref. ||13||). 



However, in general, B has been considered as an effective parameter, and 
the bag model should be only regarded as a simple model with limited con- 
nection to QCD. [] We will thus use the bag model, but only to compare our 

1 The point is that QCD sum rules [^l) lead to a bag constant, which is substantially 
larger than the value fitted in the simple bag model. This theme is stressed, e.g., in ref. || 
pag. 367: "There do exist a number of modified bag models with a smaller bag radius and 
correspondingly larger B values . . . and the authors claim that the sum rule values would 
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results with it qualitatively. 

If the compact objects in Her X-l and 4U 1820-30 are really strange stars, 
then there will be very deep consequences for both the physics of strong 
interactions and astrophysics. The existence of stable nugget of SM would 
have also profound implications in cosmology, in relation to the dark matter 
problem, in the sense that many such stars, unlike the two mentioned above, 
may as yet be unobserved. Note however, that by strange stars now we mean 
conglomeration of chargeless beta-stable (u,d,s) quarks held together strongly 
by gravity, with a surface density about five times that of nuclear matter and 
an interior density which may be 3 times higher. Witten's original conjecture 
is still applicable for such conglomerates, - that they may remain as relics 
of the cooling of the universe a millionth of a second after its formation [|]]. 
This is largely because chargeless beta-stable (u,d,s) quarks have a minimum 
in its EOS at about five times that of nuclear matter and can thus form the 
surface of a SS. 

Motivated by the fundamental importance of this issue, and by the criti- 
cism to the bag model mentioned above, in the present work we investigate 
the properties of SS, using an EOS for SM based on an alternative to the 
bag model. 

Intuitively we know that at high density the quarks should go over from 
their constituent masses to their current masses, thus restoring the approxi- 
mate chiral symmetry of QCD. On the other hand, the interquark interaction 
should be screened in the medium. The latter will give rise to deconfinement 
at high density. To this end we use the following Hamiltonian: Q 

H = £ (a, • Pi + AM,) + £ W V(j (1) 

i i<j 

where we have two potentials: a scalar and a vector. The scalar originates 
from the mass term. The quark mass, Mj, is taken to be density dependent 

support their model. Such claims are, however, not very sensible. In each of these models 
additional interactions are introduced. The energy and pressure associated with these 
additional interactions and that due to the QCD vacuum state are interwoven, such that 
it is not clear which quantity should be compared with the trace of the energy momentum 
tensor." 

2 In the following we use units where h = c = 1 
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and of the form: 

Mi = m, + (310 MeV)sech(y— ), i = u,d,s. (2) 

Po 

where ps = [p u + Pd + Ps)/3 is the baryon number density, p = 0.17 fm~ 3 
is the normal nuclear matter density, and v is a parameter. 

At high p# the quark mass Mj falls from its constituent value to its current 
one which we take to be |IT |: m u = 4 MeV, = 7 MeV, m s = 150 MeV. 



The density dependence introduces a density dependent scalar potential and 
more importantly restores chiral symmetry smoothly at high density (Fig.l). 

It is hard to justify the density dependent mass term except qualitatively. 
It appears we are the first to suggest such a form. On the other hand at 
finite T the restoration of chiral symmetry can be established from lattice 
calculations (see for example [15[]). However, from calculations of meson 



1(| or baryon properties [17|, 13], it is clear that inside a hadron the quark 



mass varies with the radius even as the density varies with the radius. For 



example in [fLq| , the minimum quark masses are, m u = 186 MeV and m s 
364 MeV whereas the averaged effective quark masses are m u = 305 MeV 
and m s = 465 MeV. It will become clear that in our calculation we are 
not dealing with small variations in the quark density as in a bound-state 
hadron calculation, and thus this density dependence of the scalar potential 
is all the more crucial for us. 

The colour dependent vector potential in eq. ([!]) is an interquark potential 
originating from gluon exchanges, and the A-s are the color SU(3) matrices 
for the two interacting quarks. In the absence of an accurate evaluation 
of the potential (e.g. from large N c planar diagrams) we borrow it from 
meson phenomenology, namely the Richardson potential |Tj|. It incorporates 
the two concepts of asymptotic freedom and linear quark confinement. The 
potential reproduces heavy meson spectra and has been used in sophisticated 
calculations involving two-body-Dirac equation derived from Dirac constraint 
mechanism for light as well as heavy mesons fl6|]. It has been well tested for 
baryons in Fock calculations [[17], The potential used for the meson and 
baryon is |19|] 

V Un 1 1 (3) 

13 ~ "27"/n(l + (k i -k J ) 2 /A 2 ) (k t -k,) 2 ' ( } 



with the scale parameter |IIJ A = 100 MeV . This bare potential in a medium 
will be screened due to pair creation and infrared divergence. The inverse 
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screening length, D , to the lowest order is ||20|| : 

P- 1 ) 2 -^ E H^iHY + rn* (4) 

i=u,d,s, 

where k{ , the Fermi momentum of the i-th quark is obtained from the cor- 
responding number density: 

k{ = (p^r 3 (5) 

and a is the perturbative quark gluon coupling. To simplify numerical 
calculations, instead of summing over all the individual flavours, we have 
averaged over the flavours so that 



where 



(D- 1 ) 2 * 6 -^^kfJ(kLr + ml v (6) 

7T 

^ = ( ,w /3 = ( M3+(fcfl3+M3 ) 1/3 . 



i) 3 

3 

m u + m d + m s 



(7) 

mav = ~ • (8) 

The model we adopt has the following attractive features: It truly describes 
deconfmed quarks at finite density, through the Debye screening (DS, in 
short). At zero density (for an isolated hadron) 

D~ x vanishes, leading to confinement. At finite density, due to DS, the 
gluon polarization acquires a non-zero value leading to deconfinement. The 
scalar potential also decreases with density, thus restoring chiral symmetry 
at high density. The resulting inverse DS is also plotted in Fig.l. 

Within the present approach, the energy density of SM can be written 

e = s k + s v (9) 
where the kinetic part is given by 



i=u,d,s % 

(10) 
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and the potential contribution is given by 

£ v = -7TT E / +1 rfx /" J fe / t k i 2 f(K kj, Mi, Mj, x)V[D-\ (ki - k,-) 2 ]dM*i 

Z7r° 77 7-1 JO JO 

(11) 

where V is the screened Richardson potential. In other words, (kj — k,) in 
eqn. @ is replaced by [(kj — kj) 2 + D~ 2 } and 

k 2 k 2 1 

f{h, kj, M u Mj, x) = ( ei . ej + 2.k i .k j .x + — - — (12) 



(e i -M i )(e i -M i ) 



with 



l k 2 i +M 2 + M l (13) 

Then we calculate the self consistent chemical potentials to satisfy the j3- 
equilibrium and charge neutrality conditions, 

= , /J-d = + He ; (14) 

2(*J) 8 - (H) 3 - {k{f - {k{f = (15) 

where /i-s are the chemical potentials of u,d,s quarks and the electron, e. We 
assume that the neutrinos have left the system (//„ = 0). 
If m e is the electron mass, k{ is obtained through 



k{ = - m\ . (16) 

To achieve the conditions (0,|T^), we have the additional difficulties origi- 
nating from the density dependence of the quark masses and DS length. 
The chemical potential for the i-th quark is given by : 



Hi= v / M2 + (fc/) 2 + (A/i i ) M + (A/i i )y, i = u,d,s. (17) 

where (A/i^m is the contribution from Ek (eq.[L(]), which is evaluated straight- 
forwardly: 

(An ) _ 9e k dM, _ d dM, 

(A/ii)M = dW^p~ ~ {ps +Ps+ Ps) ~dp~ (18) 

p\ being the scalar density for the i-th quark 



(k{) 2 + M 2 + k{ 



2M t k f J(k{) 2 + M 2 - 2Mfln ^ y U — * «■] (19) 
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On the other hand, the contribution from the potential part (|TT| ) of the 
energy density 

de v dk{ 3e v dM { 

is rather complicated and it is evaluated numerically. Equations ( HjTjp can 
now be satisfied and the EOS for /3-stable SM is obtained. 

With this EOS we solve the Tolman-Oppenheimer-Volkov equation to cal- 
culate the structure of non-rotating SS. The properties of the maximum mass 
configuration for different choices of our model parameters are summarized 
in Tab.l. Our EOS is most sensitive to the parameters v and «o which rule 
the density dependence of the quark mass and DS length respectively. How- 
ever, we found that a change within 20% of these parameters, e.g., around 
the values v = 0.33, a = 0.20, produce a change of M max and of the corre- 
sponding radius which is smaller than 10%. In the same table we also report, 
for comparison, the results we obtained using an EOS for SM based on the 
bag model with the following choice of the parameters: B = 110 MeV j fm 3 , 
m s = 150 MeV, and B = 110 MeV/ fm?, m s = 0, with the strong coupling 
constant a c = (no gluons) in both cases. There have been other theoretical 



approaches to the study of SS P, 21 1. The work of Drago et al. M uses 



(u,d) quark masses close to 100 MeV in the colour dielectric model and gets 
a maximum mass ~ 1.6 M/M Q and a radius of about 10 km. In [2I| a 



very strong magnetic field is proposed as a mechanism for softening a bag 
equation of state. 

The calculated mass-radius relations are plotted in Fig. 2, and refer to 
our present model (curves labeled eosl and eos2) and to the bag model for 
the EOS. There is qualitative agreement with the bag model results, and in 
particular for low values of the mass, M is proportional to R 3 . However, 
the two bag model calculations, close to our curves indeed use a large bag 
pressure which remains constant at all densities. In Fig. 2, we also compare 
the theoretical MR relations with the semiempirical MR for the two strange 
star candidates. The closed region in Fig.2, labeled 4U 1820-30 represents 
the semiempirical MR relation recently extracted from observational data 
by Haberl and Titarchuk [1C] and used in the theoretical analysis of ref.[Q]. 
The trapezium-like region labeled Her X-l represents the semiempirical MR 
relation for the compact object in Her X-l. We followed the analysis by 
Wasserman and Shapiro p|J updated with new mass and distance measures 



7 



of Her X-l as reported in ref. |23], 2l|. Dashed curves a and b in Fig. 2 denote 
the MR relation for Her X-l, assuming an x-ray luminosity L = 2.0 • 10 37 
erg/s and L = 5.0 • 10 37 erg/s respectively. The two values of the luminosity 
we used correspond to a distance of 5.0 kpc and 7.9 kpc respectively. To 
be aware of the sensitivity to the luminosity, we also report the MR relation 
(curve c) assuming the Eddington luminosity for spherical accretion as an 
upper bound. The two horizontal lines locate the measured mass of Her X-l 
(1.1 - 1.8) M . 

It is very important to stress that the above MR relations have been 
extracted from two different type of astronomical phenomena — x-ray burst 
spectra (4U 1820-30) and cyclotron line data from a x-ray pulsar in a binary 
system (Her X-l) — and using different theoretical models to analyze the 
original observational data. The two semiempirical MR relations overlap in 
a region of the MR plane indicating the existence of a compact object with 
a radius of 6-8 km. This shows that the analysis performed in ref. P] in 
the case of the x-ray buster 4U 1820-30 also extends to the case of Her X-l. 
In particular, neutron star models based on "conventional" EOS of dense 
matter are unable to reproduce the semiempirical MR relation for these 
two compact objects. 

The MR relation calculated with our EOS for SM is well within both 
the semiempirical MR relations of the two SS candidates. In conclusion, our 
present results confirm and reinforce the recent claimjl], that the compact 
objects associated with the x-ray pulsar Her X-l, and with the x-ray burster 
4U 1820-30 are strange stars. 

This paper was inspired by discussions of an early work done by two 
of us (JD and MD) with Prof. Jean Le Tourneux in 1984 (p7|). JD and 
MD acknowledge the hospitality of the International Centre of Theoretical 
Physics, Trieste, Italy where this work was inititated. MD is particularly 
grateful to Dr. M. Malheiro for the introduction of the subject through 
series of e-mails, specially, the units and the TOV equation. It is a pleasure 
to thank Drs. S. K. Samaddar, A. Ganguly, B. Datta and F. Weber, S.Raha 
and S. Ghosh for many discussions. 
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Table 1: Properties of the maximum mass strange star configuration obtained 
for different equations of state: Mq is the gravitational (maximum) mass, 
R is the corresponding radius, p c the central number density, e c the central 
mass density. Our EOS for different choices of the parameters are denoted as 
follow: (eosl) v = 0.40, a = 0.20; (eos2) v = 1/3, a = 0.20; (eos3) i/ = 1/3, 
cto = 0.15. The the bag model EOS are denoted as follow: (B110;150) 
B = 110 MeV/fm 3 , m s = 150 MeV; (B110;0) B = 110 MeV/fm 3 , m s = 
(a c = in both cases) 



EOS 


Mq/Mq 


R(km) 


Pc(fm 3 ) 


e c (10 u g/cm 3 ) 


eosl 


1.414 


7.18 


2.24 


AA.Q 


eos2 


1.286 


6.53 


2.57 


53.6 


eos3 


1.301 


6.62 


2.55 


52.4 


B110;150 


1.358 


7.50 


1.87 


42.0 


B110;0 


1.448 


7.90 


1.76 


38.1 



Fig. 1 The inverse Debye screening length D^ 1 for a = 0.20, and the d 
quark mass M d , for v = 0.4 (dashed line) and u — 1/3 (continuous line), are 
plotted as a function of the baryon number density. 

Fig. 2 The theoretical MR relations (curves eosl and eos2) calculated 
within the present model are compared with the semiempirical MR rela- 
tions for 4U 1820-30 and Her X-l (closed regions in the MR plane labeled 
4U 1820-30 and Her X-l, respectively). The remaining two continuous curves 
represent the MR relation calculated with the bag model EOS, with m s = 
(upper curve) and m s = 150 MeV (lower curve), and B = 110 MeV/ fm 3 in 
both cases. 
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